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Abstract. 

We present a model for inclusive charged lepton-nucleon and (anti)neutrino-nucleon cross sec- 
tions at momentum transfer squared, Q^, ^ 1 GeV^. We quantify the impact of existing low-Q 
charged-lepton deep-inelastic scattering (DIS) data on effects due to high-twist operators and on the 
extraction of parton distribution functions (PDFs). No evidence is found for twist-6 contributions to 
structure functions (SF), and for a twist-4 term in the logitudinal SF at x > 0.1. We find that DIS 
data are consistent with the NNLO QCD approximation with the target mass and phenomenological 
high twist corrections. For < 1 GeV^, we extend extrapolation of the operator product expan- 
sion, preserving the low-Q current-conservation theorems. The procedure yields a good description 
of data down to ~ 0.5 GeV^. An updated set of PDFs with reduced uncertainty and applicable 
down to small momentum transfers in the lepton-nucleon scattering is obtained. 
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INTRODUCTION 

At high momentum transfer Q the lepton-nucleon cross-sections are well described in 
terms of parton distributions (PDFs), whose evolution is well-understood in pertur- 
bative Quantum Chromo-Dynamics (QCD). The universality of the partonic description 
allows to obtain predictions for a variety of probes (e, jU, v,v) and targets, which have 
been extensively verified by experiments. However, by lowering progressively Q non- 
perturbative phenomena become more and more important for a precise modeling of 
cross-sections. 

The existing DIS data at small momentum transfer Q in principle can shed light on the 
interplay between perturbative and non-perturbative phenomena and clarify the limits 
of applicability of the parton model. Furthermore they can put valuable constraints on 
the parton distribution functions (PDFs) and on the strong coupling constant Us, due 
to their excellent statistical precision. However, significant high-order QCD corrections 
are required to study such kinematical region. The recent progress in the NNLO QCD 
calculations [1] allows to use the DIS data down to 2 ~ 1 GeV in global QCD fits by 
keeping the perturbative stability under control. 



In the formalism of the Operator Product Expansion (OPE) unpolarized structure 
functions can be expressed in terms of powers of 1/2^ (power corrections): 

F2.r,(^.Q')-F^;l(x.Q^) + "2^ + "^^ (1) 

The first term (t = 2), expressed in terms of PDFs, represents the Leading Twist (HT) 
describing the scattering off a free quark and is responsible for the scaling of SF via 
perturbative QCD (Xs{Q^) corrections. The Higher Twist (HT) terms (t = 4,6) reflect 
instead the strength of multi-parton correlations (qq and qg). Since such corrections 
spoil factorization one has to consider their impact on the PDFs extracted in the analysis 
of low-2 data. Due to their non-perturbative origin, current models can only provide 
a qualitative description for such contributions, which must then be determined phe- 
nomenologically from data. 

Existing information about high twist terms in lepton-nucleon structure functions 
is scarce and somewhat controversial. Early analyses [0, 0] suggested a significant 
HT contribution to the longitudinal SF Fl. The subsequent studies with both charged 
leptons [0, [si 0] and neutrinos [7] raised the question of a possible dependence on the 
order of QCD calculation used for the leding twist. The common wisdom is generally 
that HTs only affect the region of 2^ ~ 1 3 GeV^ and can be neglected in the extraction 
of the leading twist. 

In this communication we report our results on using the DIS data down to 2 = 1 GeV 
in the global QCD fit of PDFs with power corrections included in the analysis. 



PROCEDURE 

The analyzed data set consist of the world charged-leptons DIS cross section data for 
the proton and deuteron targets by the SLAC, BCDMS, NMC, FNAL-E-665, HI, and 
ZEUS experiments supplemented by the fixed-target Drell-Yan data, the latter constrain 
the sea quark distribution, which is poorly determined by the DIS data alone. Basically 
the same combination of data was used in the earlier fit of Ref.[|9|] with the cut > 
2.5 GeV^ imposed on the DIS data. In the present fit alongside with the softer cut 
imposed on the SLAC and NMC data, Q > 1 GeV, we also add the DIS data by 
FNAL-E-665 experiment [[lo|l since they give additional constraint on the PDFs at 
small X provided not too stringent cut on Q is applied. The cut on invariant mass of 
the hadron system W > l.S GeV is imposed on the DIS data to avoid the resonance 
region. The total number of data points (NDP) used in the fit is 3076, in the range of 
X = 0.0001 ^0.9. The analysis Q is performed in the NNLO QCD approximation with 
the target mass corrections ^ taken into account and the dynamical twist-4 (twist-6) 
terms parameterized in the additive form as model independent spline functions H{x). 
Deuteron data are corrected for nuclear effects following the model [ TH ] . 



Details of the theoretical ansatz can be found in Ref.f'o'l. 
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FIGURE 1. Left figure: the 1(7 error bands for the twist-4 (solid lines) and twist-6 (dashes) terms in 
the proton F2 (upper panel) and Fj (lower panel). The arrow indicates the region of x with the limited 
potential for the determination of twist-6 terms due to the cut on W. Right figure: pulls corresponding 
to the fits with and without twist-6 terms. The arrows in the upper panel indicate the high-y and low-j 
regions for the SLAC and BCDMS data. The data points for the JLAB-E-118 experiment at g 1 GeV 
were not used in the fit. 



In addition, the recent neutrino and antineutrino cross-section data from the CHORUS 
experiment [12] are added to the global fit for Q > l.O GeV and x > 0.045, mainly to 
constrain the corresponding HT terms. 

DETERMINATION OF HIGH TWIST TERMS 

As a first step of the analysis we checked possible twist-6 contributions to the DIS 
structure functions by keeping the terms H'''^^{x) in Eq.([Tl). Due to the W cut, twist- 
6 terms are insensitive to the large-x data and therefore we set them to zero at .x: > 0.5. 
Values of the HT terms atx — were also set to zero in view of the fact no clear evidence 
of saturation effects was found in the small x HERA data. 

Figure \T\ (left) shows the HT terms in F2 j obtained in such a variant of our fit. 
Surprisingly, we observe a large positive contribution from the twist-6 term to Fj at 
X ~ 0. 15. At the same time this contribution is compensated by a corresponding negative 
twist-4 term. Since the twist-4 and twist-6 coefficients are similar in magnitude but 
opposite in sign, the overall sum of HT terms demonstrate a weak dependence on Q. 
This observation leads us to the conclusion that in fact the twist-6 term in Fj absorbs 
some non-power- like effects. 

When we remove twist-6 terms from our fits the magnitude of the pulls appears to be 
maximal not at the lowest Q values, but at 2^ ~ 7 GeV^, exactly in the region of overlap 
between SLAC and BCDMS data (see Fig {!] right). Indeed, these two data sets show a 





FIGURE 2. Left figure: the 1 a error bands for the high-twist terms in the isospin-symmetric combina- 
tions of different structure functions (soHd lines: F2, dashes: Fj, dots: Fi) for charged leptons. Right figure: 
corresponding 1 a bands for neutrino scattering off an isoscalar target (upper panel: F2, lower panel: xFj). 
The predictions for F2 from charged leptons rescaled by the corresponding leading twist terms are also 
shown for comparison. 



certain discrepancy, which is directly translated into a fake twist-6 contribution, if such 
term is fitted as well. The values of inelasticity y are large for SLAC data at the highest Q 
and small for the BCDMS data at the lowest Q. For this reason, the discrepancy affects 
mainly Fj, which is more sensitive to y than F2. On the other hand the SLAC/BCDMS 
inconsistency at ^ 1 GeV^ is not crucial for the determination of twist-4 terms. If 
we rescale the uncertainties of data in this region to bring the pulls at the level of Icr, 
we observe a negligible increase in the corresponding HT errors. Evidently, HT terms 
are driven by data at the lowest Q available, so that in principle the SLAC/BCDMS 
data around ^ 1 GeV^ can even be dropped without any loss of statistical power. 
Preliminary data from the experiment JLAB-E-118 agree with the SLAC data at 
low Q, thus favouring the reliability of the latter at 2 ~ 1 GeV, regardless of potential 
problems in the region of overlap with BCDMS. 

Coming to the conclusion that the twist-6 terms observed are just an artefact due 
to certain inconsistencies in the data, we drop them from the final version of the fit. 
This results in a value of ;^^/NDP of 3815/3076=1.25, which is higher than unity in 
consideration of the data discrepancies discussed above. Nonetheless the problematic 
data points are spread out more or less randomly over kinematics and they do not bias 
the results of the fit. A rescaling of the uncertainties in the data points with the largest 
pulls, such that the overall ;f^/NDP becomes unity, leads to a modest increase in the 
errors of PDFs and HTs within 20%. Figure [21 (left) shows the twist-4 terms obtained 
in the final version of our fit. The HT contribution to Fj is remarkably similar to the 
one in F2, despite the two terms were fitted independently. As a result the HT term in 
Fi defined as //l = ^2 —Hj is well comparable to zero within the uncertainties. In the 



final version of the fit we then impose the constraint H2 = Hj for the isospin-symmetric 
combinations of structure functions. 

Our results indicate the HT contribution to the structure function R = Oi^/ot is also 
small in the whole range of x considered. This is in contrast to the conclusion of Ref.[0] 
about a large HT contribution to R. We explain such difference by the fact that the low-Q 
part of the SLAC data was not considered in Ref.[2]. The extrapolation of those results 
to the lower Q value must be in disagreement with data. In Fig. [3] (left) we compare our 
predictions for R with the empirical parameterization i?sLAC from SLAC data yO. The 
latter is consistent with our calculation based upon the fit with fake twist-6 terms. This 
indicates the i?sLAC parameterization is the result of the same inconsistency between 
SLAC and BCDMS data we discussed above. 

It is interesting to determine the HT contributions to (anti)neutrino structure func- 
tions independently from the ones extracted from charged lepton data. Due to the struc- 
ture of the weak Charged Current (CC) some similarities with charged leptons could 
be expected for F2 and Ft. Since the target nucleon is mostly isoscalar for neutrino 
data we impose the constraint H2 = Hj, consistently with charged leptons E|. Figured 
(right) summarizes our observations. The ratio i^^J^/Fj^^ is remarkably similar for both 
(anti)neutrinos and charged leptons over the entire x range. In addition, the use of neu- 
trino data allows us to extract simultaneously the HT contribution to xFj,. Our results 
indicate overall provides a negative contribution to the Gross-Llewellyn-Smith inte- 
gral, which is consistent with the predictions of Ref. [|l4|] . 

The total contribution of the HT terms into the DIS cross section turns out to be small 
compared to the leading-twist (LT) part. For a realistic DIS kinematics the ratio of the 
HT and LT terms is < 10%, which justifies the use of the twist expansion in our analysis. 

Finally, our results demonstrate high twist contributions to unpolarized structure func- 
tions do not vanish in the NNLO approximation. Indeed we find no strong dependence 
upon the order of QCD calculation used in the leading twist. 

IMPACT ON LEADING TWIST 

The use of additional data at low Q values allows to achieve a better separation between 
leading and higher twists by exploiting their different dependence. The correlation 
coefficients are indeed substantially reduced in the whole x range by extending the lower 
cutoff from 2.5 GeV^ to 1.0 GeV^ and do not exceed 0.3. This in turn results in reduced 
theoretical uncertainties. 

Our value of a^{Mz) = 0. 1 128 ±0.001 1 is consistent with the one obtained in Ref.[@] 
with the cut > 2.5 GeV^ and is in good agreement with the result from the non-singlet 
DIS data analysis [fisl] . The PDFs obtained in the fit with the low-Q DIS data included 
are also close to the ones of Ref.[|9|]. This manifests a good separation of the LT and 



This relation does not hold in general since the presence of an axial- vector current introduces significant 
HT contributions to in the limit of small x values and vanishing Q (PCAC). This will be discussed in 
the following. However, we are here focused on a kinematic region in which such effects are marginal 
(x > 0.045 and e > 1.0 GeV). 




FIGURE 3. Left figure: comparison of the calculations for R = Gl/Gt including high twist contribu- 
tions with the empirical parameterization of SLAC data /?slac [3]. The data points represent SLAC data. 
Right figure; the 1 a error bands for the gluon distribution obtained in our fit (solid lines) compared to one 
ofRef.ltl(dashes). 



HT terms in the fit and a stability of the perturbative calculation. The most significant 
change is observed in the gluon distribution at jc ~ 0.3, which is enhanced in the low- 
Q fit (see Figl3] right). This is due to the significant twist-4 term appearing in R in the 
fit with > 2.5 GeV^, similarly to the analysis of Ref.[2]. The LT terms in R is then 
correspondingly suppressed. Since the LT in R is proportional to the value of the gluon 
distribution, the latter is also suppressed as a result. 

The uncertainties on PDFs are improved as compared to the fit of Ref.[@]. In particular, 
the J-quark distribution is now determined within few per cent at x 0.2, which is 
comparable to the precision of the a-quark distribution. This improvement has important 
phenomenological consequences for the extraction of the weak mixing angle from 
neutrino data. For instance the analysis by the NuTeV collaboration [|l6ll . based on the 
Paschos-Wolfenstein relation, requires a good knowledge of the valence distributions 
in order to guarantee an accurate correction for the target non-isoscalarity [|l7ll . This 
correction is proportional to the ratio xi/xq, where xi and xq are the integrals over x of 
the iso-vector and iso-scalar combinations of the valence quarks, respectively. From our 
fit we obtain xi/xo = 0.424(6) , which provides an uncertainty on the non-isoscalarity of 
the target comparable to the experimental uncertainties from NuTeV data. 



EXTRAPOLATION TO Q < 1 GEV 

In order to describe the structure functions in the region of low momentum transfer we 
apply a smooth interpolation between the high regime, which is described in QCD 
in terms of LT and HT contributions as discussed above, and the ^ predictions 
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FIGURE 4. Left figure: interpolation of structure functions in the region 

< < 1 GeV^. The example 

given in the plot refers to F2 for charged lepton scattering on protons at x=0.07 (see text for details). Right 
figure: pulls with respect to the calculations from our fit for charged leptons as a function of Q^. The 
region <l GeV^ and data from JLab E-118 were not included in the fit. 

derived from current conservation arguments. We use Q^^ = 1 GeV^ as the matching 
point between high- and low-Q^ regimes. 

The conservation of the vector current (CVC) suggests that for the electromagnetic 
interaction in the limit of 2^ ^ the structure function Fj vanishes as while Fi 
vanishes faster than Fj in such a way that Fl/Fj ~ Q^. Although these asymptotic 
conditions define the value of structure functions at = 0, they do not specify at which 
scale such behaviour sets in. On the other hand, at Q^^ we know precisely the value of 
SFs from the twist expansion. In the region < < Qf„ we interpolate by using cubic 
splines calculated for fixed x values. The corresponding coefficients are fully determined 
by the condition both functions and derivatives should match with the twist expansion at 
Ol^. Figure in (right) illustrates the interpolation procedure for F2 on protons in charged 
lepton scattering. 

We obtain a good description of charged lepton data down to ~ 0.5 GeV^, as can 
be seen from Fig. |4] (left) showing pulls as a function of at a fixed x value. Below 
that some deviations are visible in the recent JLab E-118 data. This may indicate the 
transition from dynamics to the behaviour dictated by gauge invariance is slower than in 
our simple extrapolation of the twist expansion. 

In the low-2 region (anti)neutrino cross sections are dominated by the longitudinal 
structure function Fi and the latter is driven by the axial-current interactions. Similarly 
to the charged lepton case, the structure function Ft vanishes as at low Q^. This 
behaviour holds for both the vector and the axial- vector contributions. However, in 
the longitudinal channel the low- 2 behavior of the vector and axial- vector parts are 
different. The vector component F^ still vanishes as Q'^ at low values. 




FIGURE 5. Left figure: the PCAC contribution to the neutrino structure function F2 calculated for 
X = 10^^ as a function of for a few different targets (labels on the curves). A value MpcAc = 0.8 
GeV is assumed. Right figure: comparison of the ratio R = Fi/Ft calculated for the isoscalar nucleon for 
the charged lepton (dashed-dotted line) and neutrino (solid line) cases at x = 0.125. Also shown are the 
results for different nuclear targets ('^C, ^^Fe and -'^^Pb from top to bottom). A value MpcAc = 0.8 GeV 
is assumed. Determinations from SLAG yj, CCFR [19,1 and CHORUS L1Z1 are given for comparison. 



In contrast to the vector current, the axial-vector current is not conserved. For low 
momentum transfer the divergence of the axial-vector current is proportional to the 
pion field (Partially Conserved Axial Current or PCAC): dA^ = fjtin^cp'^ where 
is the pion mass, fj^ = 0.93mj^ the pion decay constant and (p^ the pion field in the 
corresponding charge state. We introduce explicitely a PCAC contribution to F^: 

Ft = ^Fl^^^f^CKdQ')+Ft (2) 

where 7= (1 +Ax^M'^ /Q^yl^, fI^^^ = f^Ojt/7t and On = Onis^Q^) is the total cross 
section for the scattering of a virtual pion with the center-of-mass energy squared 
s. The last term is similar to F^ and vanishes as Q^. Since the PCAC contri- 
bution is expected to vanish at high we introduce a form factor /pcAc(2^) = 
(1 + t^^]' where the dipole form is motivated by meson dominance ar- 

guments H. It is important to note the pion pole does not directly contribute to structure 
functions and hence the mass scale controlling the PCAC mechanism, MpcAC^ cannot be 
the pion mass itself, but is rather related to higher mass states like ai, pn etc. This scale 
is not well known from theory and must be determined with data themselves. A value 
M^CAC = 0.8 GeV seems to provide the best agreement with existing low-Q data. 



Formally the PCAC contribution can be considered as a high twist term. 
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FIGURE 6. Pulls for neutrino (left) and antineutrino (right) cross-section data from CHORUS (l^ with 
respect to our calculations as a function of Q^. Nuclear corrections ifTHl for the lead target and electroweak 
corrections ll2Tll are taken into account. 



Since F2 = {Fi+FT)/Ythe presence of the PC AC terms implies the structure function 
F2 at low is dominated by the nonvanishing fI'^^^ term. Figure [5] (left) illustrates the 
magnitude of such contribution to F2 for an isoscalar nucleon and for a few nuclear 
targets [18]. The values of F2 for heavy nuclear targets are systematically smaller 
because of nuclear shadowing effect for the pion cross section. Our predictions for the 
asymptotic value are consistent with the determination F2{Q^ = 0) = 0.210 ±0.002 by 
the CCFR experiment on an iron target ll20(l . 

It is instructive to compare the low-Q^ behaviour of R = Fl/Ft for charged-lepton 
and neutrino scattering. In both cases Fj oc as 0. However, if F^ oc for the 

electromagnetic current, for the weak current Fl — > pJCAC j^j^^^ ^^^^ vanish 
in the low-Q^ limit. Then the behavior oi R dX GeV^ is substantially different 

for charged-lepton and neutrino scattering. In Figure [5] (right) we illustrate this effect 
by calculating 7? as a function of ^ at a fixed x for an isoscalar nucleon and a number 
of nuclei. We observe the nuclear correction partially reduces the differences between 
charged leptons and neutrinos and smoothes out the divergence of R in the latter case. 

Figure [6] shows the pulls of our fit to neutrino and antineutrino inclusive inelastic 
cross-section data from the CHORUS experiment fv^ on a lead target. The calculation 
includes nuclear corrections [|lll] for the lead target, electroweak corrections uM and 
the PC AC contribution. 

The determination of LT and HT terms is performed from all available data with > 
1.0 GeV^ and W > l.S GeV. It is interesting to check the extrapolation of DIS structure 
functions into the resonance region for < 1.8 GeV. Our results are consistent with the 



duality principle and the integral of the difference between the recent JLab resonance 
data and the average DIS predictions is consistent with zero within few percent in the 
entire kinematic region. This can be also considered as an indirect indication in favour 
of a negligible twist-6 contribution to structure functions. 
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